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Abstract 

We develop a new fluid model of a warm plasma that includes the radiative self-force on each plasma 
electron. Our approach is a natural generalization of established methods for generating fluid models without 
radiation reaction. The equilibrium of a magnetized plasma is analysed, and it is shown that the thermal 
motion is confined to the magnetic field lines. A dispersion relation is deduced for electric waves in a 
magnetized plasma, and it is shown to agree with our recently established relativistic kinetic theory derived 
from the Lorentz-Abraham-Dirac equation. 

1 Introduction 

Contemporary advances in ultra-intense laser facilities have driven the recent surge of interest in the collective 
behaviour of charged matter in extreme conditions, and a particularly vexing topic in that context concerns the 
coupling of an electron to its own radiation field. An accelerating electron emits electromagnetic radiation, and 
the energy and momentum carried away by the electromagnetic field must be accounted for. In most practical 
cases, the Lorentz force on an electron, due to an applied electromagnetic field, is considerably larger than the 
force due to the electron's emission and the effect of the recoil due to the emitted radiation is negligible or 
can be adequately represented using simple physical reasoning. Although such arguments avoid the difficulties 
that plague more comprehensive analyses, the parameter regimes promised by forthcoming ultra-intense laser 
facilities ensure that more fundamental considerations are now of practical necessity. For example, ELI |T] is 
expected to operate with intensities lO^^W/cm^ and electron energies in the GeV range, at which level the 
radiation reaction force becomes comparable to, and can even exceed, the applied force due to the laser field. 

One of the most notorious equations in physics was developed by Dirac f2] to describe a classical point 
electron's radiative self-force. The Lorentz-Abraham-Dirac (LAD) equation is a fully relativistic description of 
a structureless point particle in an applied electromagnetic field Fab and has the form 

with q the particle's charge, m the particle's rest mass, r = q^ /Q-Kra in Heaviside-Lorentz units with c = = 
Ho = 1, and the tensor is 

For an electron, q = — e < 0. The Einstein summation convention is used throughout the present article, indices 
are raised and lowered using the metric tensor rjab = diag(— 1, 1, 1, 1) and lowercase Latin indices range over 
0, 1, 2, 3. The particle's 4-velocity dx°'/dX is normalized as follows: 

dx dxa . , 
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where A is the particle's proper time. 

Dirac derived ([l} for an electron by appealing to the conservation condition on the stress-energy-momentum 
tensor (see Ref. [3] for a recent discussion of the derivation). Dirac's approach requires a regularization of 
the electron's singular contribution to the stress-energy-momentum tensor followed by a renormalization of the 
electron's rest mass. The procedure leads to the third-order term in ([1]), which is the source of the famous 
pathological behaviour exhibited by solutions to the LAD equation (see Ref. [3] for a recent discussion). 

The standard approach to ameliorating the problems with the LAD equation is to replace the third-order 
terms in ([T]) (radiation reaction force) with the derivative of the first term on the right-hand side of ^ (the 
applied Lorentz force). This procedure is justifiable if the radiation reaction force is a small perturbation to the 
Lorentz force, and it yields the Landau-Lifshitz (LL) equation: 

d^x" _ q dx'' q ^ a dx'' dx^ ,^q\a pc dx'' 

Unlike the LAD equation, the LL equation is second order in derivatives in A and its solutions are free from 
pathologies. 

Recent years have seen a substantial growth of interest in kinetic theories incorporating radiation reaction 
(see, for example, Refs. O [B]), and almost all such theories are based on the LL equation from the outset. 
However, the briefest of glances at ([T]) and ^ suggests that a degree of mathematical clarity (with concomitant 
physical insight) is likely to be gained by starting with a kinetic theory based on the LAD equation. Of course, 
procedures for extracting physically acceptable behaviour from the LAD kinetic theory must be introduced 
during the analysis. 

Remarkably, until recently [7] , very little has appeared in the literature concerning a fully relativistic many- 
body system governed by the LAD equation. In addition, some of what can be found [5] appears to lead to 
physically inconsistent results (see Ref. [7] for a discussion of this point). Perhaps one of the reasons why the 
LAD equation appears to have been neglected in the physics literature is the necessary, but unconventional, 
introduction of the notion of a 'phase' space encoding acceleration as well as velocity and spacetime events. By 
employing the appropriate geometrical machinery, we recently showed how to construct a Vlasov equation 
based on the LAD equation. 

The purpose of the present article is to demonstrate how the Vlasov equation developed in Ref. [7] naturally 
induces a warm fluid theory and to explore the consequences of that theory. Although we followed a similar 
approach in Ref. [9J, the warm fluid approximation used there did not treat the acceleration and velocity 
moments on an equal footing. Our aim in the following is to redress that balance and obtain a more general 
warm fluid theory. We also use our new theory to explore electric waves in a magnetized plasma. 

2 Kinetic theory 

As shown in Ref. [7 , the LAD equation may be written in first-order form as 

(5) 



(6) 



dx" 
H 
dv'' 
'dX 

^ = x"£av^ + - f + ^F^'axA (7) 
dX T \ m J 

in terms of the proper-velocity (celerity) v — {v^,v'^,v^) and acceleration a — (a^,a^,a'^) in a 10-dimensional 
'phase' space {x"',v^,a''') where lowercase Greek indices range over 1,2,3. The shorthand — y/l + v^, 
xP = a^Vfj,/y/l + v^, xf^ = v^^, x^ — have been used, and the parametrization in terms of v and a has been 
chosen to satisfy the constraints x'^Xa = —1, x°'Xa = arising from ^ and its first derivative with respect to 
A. Greek indices are raised and lowered using the Kronecker delta 5^^. 

From a geometrical perspective, the Vlasov equation may be understood as the preservation of a differential 
form of maximal degree along the fiow of single-particle orbits. For present purposes, the 1-particle distribution 
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/ = f{x,v,a) of plasma electrons may be naturally extracted from that differential form [7| leading to 



where L is the Liouville operator 



■ a u 9 

L = x"- ha^T^ 



Lf + -/ = 

T 



1 



q 



d 



(8) 



(9) 



The second term on the left-hand side of ([5]) may be understood as a consequence of losses due to radiation (for 
more details see Ref. [7j). 
Maxwell's equations are 



daFbc + dcFab + dbFca = 0, 

with da = d/dx"' and N"^ the number 4-current of the plasma electron fluid, 



N^ix) 



x''f{x,v,a) 



1 + 1>2 



(10) 
(11) 



(12) 



where v"^ = v^Vf^. The presence of the factor l + v"^ in ([T^ and the second term in ([5]) are related, as discussed 
in Ref. [7]. 

For clarity of presentation, the ions will be prescribed as a homogenous background and included in the 



external source 4-current Jexf 



3 Fluid theory 

Kinetic theories are not always the most convenient tools for analytical investigation of the collective dynamics 
of charged matter (whether or not they include the radiative self- force). Furthermore, extensive computational 
resources are usually required to numerically solve the integro-differential systems of equations found in such 
theories. 

In practice, the 1-particle distribution / will usually contain more information than is needed and one may 
replace the Vlasov equation with a fluid theory that encodes / using a subset of its velocity-acceleration moments. 
In general, macroscopic fluid theories are more analytically amenable and less computationally demanding than 
their kinetic counterparts, and their relationship with experiment is more immediate. 

Using the notation first introduced in Ref. H], the natural moments of / are tensor fields on spacetime 
defined as 

^ai...a,:fci...fc„(^) ^ / . . . i"' i''^ . . . i*-" / (x, -y, a) ^^^^ . (13) 



1 + 

Indices associated with velocity are located to the left of the colon in so.i---'ie-bi---bn ^ whereas those associated 
with acceleration are located to the right of the colon. An absence of indices will be denoted by as follows: 

s'-jf'^. (») 

S'^'-"'-^ = j x"^ ...x'^^f^^, (15) 

S^-^^-^- ^ j i^K..i^-f^^. (16) 

Natural moments with an immediate physical interpretation include the number 4-current A^" = 3°"''^ and 
stress-energy-momentum tensor T°-^ ~ mS""^-" of the plasma electron fluid. The scalar field is the relativistic 
enthalpy. 
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The Vlasov equation (|8]) may be recast as the infinite hierarchy of tensor equations 

(17) 
(18) 

.pb^gc-M j ^ 0, (19) 





= 0, 










daS''-^ - 


- s'--', - 






i gb:c 


- s''--'' = 0, 




_ g<l:bc 


obc-.d 
- o d-' 




cb:cd 
~ J < 


Qc:hd 
i — ^ d- 



m 



(20) 

Sb-'^ + ±F''dS'"^A ^0, (21) 
m J 

2g$:bc ^ S.f\ S'^'-" + ^F-'d S"^^^) = 0, (22) 

(23) 



TO TO 



na :0 
i-" a 




^"a = 








>-> a 


= -s'-- 


oab: 
»^ a 


= 0, 


qcl: 
O ah 


= 0, 



on spacetime, where . . . indicates equations whose first term Q^^so-i - ai-bi...b„ satisfies £ + n > 3. Furthermore, 
the identities x°'Xa = and ±"'±a = — 1 lead to the constraints 

(24) 
(25) 
(26) 
(27) 
(28) 
(29) 
(30) 

where . . . indicates equations containing natural moments with rank greater than 3. 

For practical purposes, a finite set of equations must be chosen from the infinite sequences (I17ll23p . (|24j|30l) 
and this may be achieved by introducing the bulk velocity U"", bulk acceleration A° 

and the centred moments 

X ("f -A>"{x)).^^ h"- - A'- (x)\ fix. V, a) (33) 

followed by the assumption that all centred moments of a particular rank or greater are negligible. 

The rank 1 centred moments R"-^, R^'" trivially vanish due to the definitions of 5"°-®, 5®-". Setting 

all centred moments of rank 2 or greater to zero is equivalent to demanding that the electron distribution has 
zero spread in velocity and acceleration, and this assumption collapses (|17ll23p . (|24ll30p to the LAD equation 
for U" (see Ref. [5] for details). 

A fluid modelling a collection of electrons whose distribution has a small, but non-negligible, spread about 
the bulk velocity and bulk acceleration is more subtle to construct. Following an approach analogous to the 
scheme introduced by Amendt [TD] for fluids without the radiative self-force, we introduce the scalar field 
e = VI + U'^Ua and hypothesize that R<^i -<^t-bi-b^ = C)(e'^+") with S'^ = 0{e"), U" = 0(e°), A" = C(e"). 
Thus, R"''-^ = 0(e2), R'^-'^ = 0{e^), R^'-''^ = 0{e^) and inspection reveals that the total number of independent 
components of (|171l22p matches the total number of independent components of the variables 5®, J7°, A°-^ i?"''-®, 
j^a-.b^ j^(ti:ab all O(e') tcrms in (fT7ll22l) set to zero. However, in general, it is not possible to find solutions 
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to the resulting field equations that also satisfy the constraints ()24ll29p with all 0(e'^) terms set to zero in those 
constraints. Instead, inspired by previous work on fluids without the radiative self- force [11) . we impose the 
weaker condition that (|24ll29p need only be satisfied to 0{e^), which leads to 

R''a-'^ + S'^{l + U''Ua)^Oie^), (34) 

R---^ + S^U''Aa^O{e''), (35) 

U''R%''^ =Oie^), (36) 

U^Ra " ^0{e^), (37) 

U''R''-a + AaR"""-^ - ©(e^), (38) 

VR'^-ab + AaR'^-b^Oie^). (39) 

A warm fluid model including the radiative self-force is obtained by setting to zero all terms that are 0{e^) in 
(I17ll22p . and solutions to the resulting system of PDEs are sought that satisfy p4ll39p . 

4 Equilibrium states in a magnetized plasma 

An investigation of p7ll22p . with all 0{e^) terms set to zero, leads to 

A" = 0, (40) 

F^bU^ = 0, (41) 

30:bc I 1 / n6:c , 1 



T 



JlVI:bc^_l^b:c^J_pc^j^M-M\Q^ (42) 



2m 

j^a-.bc ^ J_lJ^d:bpc^ ^ R'^-''F^d) - (44) 
2m 

when all of the fields are constant and, as expected for equilibrium states, (|i(71) shows that the bulk acceleration 
vanishes and the electric field vanishes in the fluid frame as indicated by (HI]). As expected, the above do not 
constrain the value of the magnetic field. 

The centred moments i?"''', r'^'"-^ ^^^y |-,g eliminated from (|42II44I) to give 

^{J■,{J■,7^}}-i[^,7^]=0 (45) 
where denotes the anti-commutator and [A, B\ denotes the commutator of the matrices A, B with 

7^ = (i^^^®), (46) 

:F^{F\). (47) 

An elegant method for solving exploits the algebraic properties of the electromagnetic field that arise 
because the electric field vanishes in the fluid frame. It may be shown that FF* = F*F = where the 
matrix F* — [e°'bcd F'^'^ /2) is the electromagnetic dual of J^, with Sabcd the Levi-Civita alternative symbol, and 
furthermore 

VV = V, (48) 
VV = V, (49) 

rv = rr = o, (50) 
v + r^i (51) 
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where I is the unit matrix and 

V= , . (53) 

tr(J-*2) ^ ^ 

Hence J^, J^* naturally induce a pair of orthogonal idempotent matrices V , V which can be used to decompose 
the tangent space on spacetime. The idempotent V projects vectors into the 2-dimensional subspace of the 
tangent space orthogonal to the bulk velocity C/" and the magnetic field lines in the fiuid frame. 

The stress-energy-momentum tensor of the electron fluid is T°-^ = mS""^'®, and the heat flux is chosen to 
vanish in the equilibrium state. Since 5'"''^® = S^VU^ + R"^'^ , the choice 

UaR"^-^ = (54) 

ensures that the heat flux vanishes and ([55]) is satisfied. It is then straightforward to use (^01 - [M]) to show 
that the left-hand side of the constraints ([55l [571 - IM)) vanish. 

Equations (|4n 154)) suggest the introduction of the idempotent H = {5^ — U°'Uh/\U\^), which projects 
vectors into the 3-dimensional subspace of the tangent space orthogonal to [7°, where \U\ = \/~U°-Ua- Since 
^ab:0 ^ ^6a:0 foUo^g ^7^ = 7^^ = 7^ and, furthermore, H J" = J"H = J" since the electric field vanishes in the 
fluid frame and Fab = —Fba- Hence, natural idempotents for constructing the solution to are HT^H and 
HPH, where the former is simply 7-", since HPH = 7-", and the latter projects vectors into the 1-dimensional 
subspace parallel to the magnetic field lines in the fluid frame. 

A suitable form for TZ isotropic around the magnetic fleld lines is 

Tl = ^HT'H + ^H-pn (55) 

TO TO 

where the scalar p_L (resp. p\\ ) is the pressure transverse (resp. longitudinal) to the magnetic field lines. Hence, 
using ([55]) and the above properties, it follows [F,TZ] — whereas {J", {F,TZ}} ^ in general. By inspecting 
(|45|). we conclude that — Q and p\\ can be chosen freely. It follows that the thermal motion is along the 
magnetic field lines only: 

Hpn. (56) 

TO 

The result (j56p has a straightforward physical explanation. An electron spiraling around the magnetic field 
lines will emit cyclotron radiation, whereas there is no emission due to inertial motion along the magnetic field 
lines. Hence, angular motion around the magnetic field lines will be damped by radiation friction and the 
electron's velocity will tend towards a constant vector parallel to the magnetic field lines. 

The background ion number 4-current is N^^^ — nion<5o where the constant njon is the proper number density 
of the ions, and S"'-^ = N" = iVi^„ follows fr om the Maxwell equation pT|) since the electromagnetic field is 
constant and J^^t = ^l^ton- The enthalpy 5® and bulk velocity are determined in terms of py, to, nion 
using (ini [Ml EH) and the property p[| = O(e^). The latter follows from and the assumption i?"^-" = ©(e^). 



5 Electric waves 

A small amplitude, spatially harmonic, perturbation to the equilibrium proper number density, directed along 
the magnetic field lines, induces an electric wave that propagates along the magnetic field lines. In the following 
the z = axis is chosen to lie parallel to the magnetic field lines in the rest frame of the ions. We denote 
t = and introduce the proper electron number density n = nion + Sn with Sn cx 3fi[exp(ifcz — iujt)]. Similar 
expressions are also introduced for the other field variables. 

A linearization of (|17ll22p . with all C(e^) terms set to zero, around one of the equilibria introduced in the 
previous section yields a dispersion relation w — uik) solved by 
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where the normahzed equihbrium temperature 9 = p\\/ {riioam) and the plasma frequency ojp = ■\/ n\anCp' jm. 
The sohition ((57)) was chosen because it reduces to a; = cjp in the hmit t — > and 61 — > 0. The dispersion 
relation also possesses unphysical solutions that diverge as r — > 0, which must be discarded. As required, the 
linearized solution to (|17II22I) can also be shown to satisfy the constraints ([Ml - [55)1 . 

The limit t of (|57)l yields a solution analogous to that of Clemmow and Willson's [H] relativistic 
generalization of the Bohm-Gross dispersion relation. The only difference is the numerical factor 3/4 in the 
relativistic correction to the frequency shift in (j57p . which arises because of the anisotropy of ]^°-^'-^ in equilibrium. 
The first damping term in (|57p is an old result (131 common to kinetic theories with radiation reaction. The 
second damping term depends on the temperature and agrees with recent results presented elsewhere as we 
will now show. 

The kinetic theory developed in Ref. 7 requires that the distribution / has the form 



= + v'-g{x, v)6^^^ {a - A{x, v)) (58) 

for the LL kinetic theory to be recovered, where g{x, v) is identified as the usual 1-particle distribution on 
event- velocity space and J'-'^^ is the 3-dimensional Dirac delta. In Ref. [7], the conventional relativistic Vlasov 
equation was obtained for g{x, v) in the limit r — >■ and a dispersion relation for electric waves without an 
external background field was derived in terms of the equilibrium distribution g{v). Specializing that dispersion 
relation to a spatially homogeneous and static distribution with zero transverse spread in velocity yields 

1 = ^ / ^ ^ (59) 

where 4* — uj^/l + — ku, and the equilibrium distribution has the form 'g{v^,v^,u) — h{u)S{v^)5{v'^) where 
d is the 1-dimensional Dirac delta. The judiciously chosen pre-factor \/l + v^ in (|58p ensures that the factor 
l/VT-M? appears in the measure in (j59p . 

The proper number density rijon of the ions and the normalized equilibrium temperature are related to h 
in the usual way: 

5°^\q„. = / VT+^Hu)^fi=, (60) 



I cqm 

J33:0 



^ion ^ion J vT 

and the heat flux of the electron fluid vanishes in equilibrium 



C33:W| 1 r J,, 

* loom _i_ I ^2^(^)_f^ (g^) 



du 



T'\^m = m53O^0|eq„, = m / uVlT^h{u)^== = 0. (62) 



u^ 



Hence 

riion = y h{u)du, (63) 

u^h{u)du + . . . , (64) 

u h{u)du — (65) 



where . . . indicates contributions that arise from 0{u'^) integrands. 
Using (Pll^ and 



^-2(1 + ir^')Vl + 



1 + -M^ + 5- - ZT- 1 + — + - + Oir^u-") (66) 
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it follows that (j59p yields the dispersion relation 



1 = 



2 \ u)'^ 



(67) 



where . . . indicates contributions that arise from 0{u^) terms in integrands or are ©(t^). 

Terms that arise from 0{v?) integrands are negligible in the warm fluid approximation, and it may be shown 
([57]) satisfies (p7)) when the terms indicated by . . . are dropped from the calculation. Hence, the warm fluid 
theory presented here and the results of the kinetic theory developed in Ref. [7] agree, lending support to the 
methods established herein. 



6 Conclusion 

We have introduced a natural generalization of the standard warm fluid approximation that closes the hierarchy 
of velocity-acceleration moments generated from ([S]). Our model yields a relativistic generalization of the 
Bohm-Gross dispersion relation for electric waves propagating parallel to an ambient constant magnetic field 
in a magnetized plasma. Although the calculation in Ref. [7 assumes that the background electromagnetic 
field vanishes, the results found there are applicable to the present case if the equilibrium distribution models 
electrons confined to move freely along the background magnetic field lines. The dispersion relation for electric 
waves found in Ref. agrees with the present work. 

Inspection of the imaginary part 3(a;) of u in (1571) shows that plane waves with higher k are damped less 
than those with lower k. The higher frequency components of a localized disturbance 5n in the electron number 
density n are longer-lived than the shorter frequency components. Although Sn will spread due to the k^ term 
in the real part 3?(aj) of lu in (1571) . the radiative self- force encourages additional spreading due to the relative 
enhancement of the higher frequency components. 

Further ramifications of the warm fluid theory introduced here will be presented elsewhere. 
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